Purpose: derive metrics of fragmentation and connectivity for reporting units
Load spatial data
Code
# basin
basin <- vect("/Users/jeffbaldock/Library/CloudStorage/GoogleDrive-jbaldock@uwyo.edu/Shared drives/wyo-coop-baldock/UWyoming/Snake River Cutthroat/Data/Spatial/Basin Delineation/BasinDelineation/MajorBasins_Watersheds.shp")
basin <- subset(basin, basin$site %in% c("SnakeHeadwaters"))
# reporting unit watersheds
watersheds <- vect("/Users/jeffbaldock/Library/CloudStorage/GoogleDrive-jbaldock@uwyo.edu/Shared drives/wyo-coop-baldock/UWyoming/Snake River Cutthroat/Analyses/Snake River GSI Quarto/Landscape Covariates/Watershed Delineation/RepUnits_Watersheds.shp")
# waterfalls
# falls <- terra::vect("/Users/jeffbaldock/Library/CloudStorage/GoogleDrive-jbaldock@uwyo.edu/Shared drives/wyo-coop-baldock/UWyoming/Snake River Cutthroat/Analyses/Snake River GSI Quarto/Landscape Covariates/Barriers/SARP Waterfalls/Archive_Falls_Mar082024.gdb")
# falls <- project(falls, crs(basin))
# falls <- crop(falls, basin)
# writeVector(falls, "/Users/jeffbaldock/Library/CloudStorage/GoogleDrive-jbaldock@uwyo.edu/Shared drives/wyo-coop-baldock/UWyoming/Snake River Cutthroat/Analyses/Snake River GSI Quarto/Landscape Covariates/Barriers/SARP Waterfalls/SARP_Waterfalls_SnakeHeadwaters.shp")
falls <- vect("/Users/jeffbaldock/Library/CloudStorage/GoogleDrive-jbaldock@uwyo.edu/Shared drives/wyo-coop-baldock/UWyoming/Snake River Cutthroat/Analyses/Snake River GSI Quarto/Landscape Covariates/Barriers/SARP Waterfalls/SARP_Waterfalls_SnakeHeadwaters.shp")
falls$BarrierType <- "Waterfall"
mapview(st_as_sf(falls))
Code
# barriers
barriers <- read_csv("/Users/jeffbaldock/Library/CloudStorage/GoogleDrive-jbaldock@uwyo.edu/Shared drives/wyo-coop-baldock/UWyoming/Snake River Cutthroat/Analyses/Snake River GSI Quarto/Landscape Covariates/Barriers/aquatic_barrier_ranks_Jan2025/aquatic_barrier_ranks.csv") %>% filter(Removed == "no", Passability %in% c("Complete barrier", "Partial passability", "Seasonably passable - non salmonid", "Seasonably passable - salmonid"), HasNetwork == "yes")
barriers <- vect(barriers, geom = c("lon", "lat"), crs = "+proj=longlat +datum=WGS84")
barriers <- project(barriers, crs(basin))
# combine falls and barriers
barriers_comb <- union(falls, barriers)
barriers_comb$BarrierType <- recode(barriers_comb$BarrierType, "Assessed road-related barrier" = "Road")
# map it
mapview(st_as_sf(watersheds)) + mapview(st_as_sf(barriers_comb), zcol = "BarrierType")